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We present a functional renormalization group (fRG) treatment of trigonal grapliene nanodiscs 
and composites thereof, modeled by finite-size Hubbard-like Hamiltonians with honeycomb lattice 
structure. At half filling, the noninteracting spectrum of these structures contains a certain number 
of half-filled states at the Fermi level. For the case of trigonal nanodiscs, including interactions 
between these degenerate states was argued to lead to a large ground state spin with potential 
spintronics applications [I]. Here we perform a systematic fRG fiow where the excited single- 
particle states are integrated out with a decreasing energy cutoff, yielding a renormalized low-energy 
Hamiltonian for the zero-energy states that includes effects of the excited levels. The numerical 
implementation corroborates the results obtained with a simpler Hartree-Fock treatment of the 
interaction effects within the zero-energy states only. In particular, for trigonal nanodiscs the 
degeneracy of the one-particle-states with zero-energy turns out to be very robust against infiuences 
of the higher levels. As an explanation, we give a general argument that within this fRG scheme 
the zero-energy degeneracy remains unsplit under quite general conditions and for any size of the 
trigonal nanodisc. We furthermore discuss the differences in the effective Hamiltonian and their 
ground states of single nanodiscs and composite bow-tie-shaped systems. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

Graphene-nanodiscs (GNDs) are nanostructures con- 
sisting of a finite bipartite honeycomb-lattice. Among 
them a large variety of shapes is possible. Of particular 
interest are GNDs with a large ground state degeneracy 
where interaction effects can lead to the formation of a 
high spin state with relatively long lifetime that could be 
used in spintronics applications [1- 3J . In a tight-binding- 
description metallic GNDs with half-filled zero-energy- 
states are very rare [4]. As shown in [5] the emergence 
of zero-energy-states is related to the morphology of the 
honeycomb-lattice. The number of these states 77 is equal 
to the difference 77 = a — /3, where a and /3 are the 
maximum numbers of nonadjacent vertices and edges. 
Following a classification in Ref. [B] we distinguish be- 
tween GNDs where 77 is equal to the sublattice-imbalance 
\Lb — La\ of the bipartite honeycomb-lattice consist- 
ing of the two sublattices A and B and GNDs where 
77 > \Lb — La\- La and Lb are the numbers of lattice 
sites on sublattice A and B. 

One example for the first class are trigonal zigzag- 
GNDs (cf. Fig. [ija) which are characterised by the size- 
parameter N . The sublattice-imbalance is = Lb — La 
and 77 is equal to N . In contrast, bow-tie-shaped nanos- 
tructures (cf. Fig. [l]b) represent the second class with 
zero sublattice-mismatch but with 77 > 0. 

To describe electron-electron-interactions in graphene 
nanodiscs it is common to take a p^-band Hubbard-like 



'Electronic address: |mkinza@physik.uni-wuerzburg.de| 




o sublattice A 
• sublattice B 



FIG. 1: Two different kinds of graphene-nanodiscs: (a.) Trig- 
onal zigzag-nanodiscs, that can be characterized by the size- 
parameter N . Here the zero-energy-degeneracy r; is equal to 
the sublattice-imbalance Lb — La = N (h.) bow-tie-shaped 
nanostructure with zero sublattice mismatch and 77 = 2. 

model of the form 

H = -t c| <,Cj> + UY^ n^^t<i + y ^^'<^^o^- 

(1) 

where the sum goes over all nearest neighboring sites 
The operators c^^ and Ci^^ create and annihilate 
electrons with spin a on site i. The nearest-neighbor 
hopping-amplitude t is of the order of 3eV. 

For nearest-neighbor interaction Vi = an exact theo- 
rem by Lieb [7] exists stating that the ground state spin 
of a repulsive Hubbard model on a bipartite lattice with 
sublattice site numbers La and L b and La + Lb even is 
equal to S = ^\Lb ~ La\- From this it can be expected 
that the electron-spins of trigonal zigzag-GNDs prefer 
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ferromagnetic order, while in bow-tie-shaped structures 
they adopt a total-spin-zero state. Of course, this argu- 
ment is restricted to onsite interactions, and one might 
wonder if the spin state changes for more general interac- 
tions. In an extended Hartree-Fock-approximation (8j for 
the zero-energy-states of trigonal zigzag-GNDs one finds 
a large negative exchange energy that gives rise to ferro- 
magnetic order for any small N. This remains valid in 
the case of non-local interactions. It also allowed the au- 
thor to estimate the spin excitation energies to be of the 
order of several hundred meVs. Hence the ground state 
spin seems to be rather robust. However, in this anal- 
ysis, the interaction effects are treated only within the 
subspace of zero-energy states. The question whether 
the empty or filled excited single-particle levels not in- 
cluded in Ezawa's treatment lift the ground state degen- 
eracy through virtual excitations. This splitting would 
then compete with the Hund's rule or exchange term. In 
addition, the interactions between the degenerate states 
will be altered by virtual processes through the excited 
levels. Quite generally, it would be desirable to derive an 
effective low-energy Hamiltonian by integrating out the 
higher excitation levels in a renormalization procedure 
rather than just neglecting these levels. 

In the following we will use the functional RG formal- 
ism to accomplish this task within reasonable approxima- 
tions. The fRG formalism has already proven useful in 
the study of two-dimensional Hubbard-like bulk systems, 
mainly for the search of instabilities in Fermi liquids |9|- 
112] . Furthermore, it has been applied in real space to 
many one- and zero-dimensional mesoscopic systems, giv- 
ing very good descriptions of boundary exponents in the 
density of states [13] and transport properties [14] . Here 
we show how one can use the fRG in order to derive an 
effective theory for the zero-energy-state-sector of GNDs. 
We test the method at two examples: the first one are 
trigonal nanodiscs as in Ezawa's papers. Interestingly, 
we find that within our fRG treatment at half filling, the 
single-particle levels at zero energy do not split up at all 
while the interaction parameters get renormalized by in- 
tegrating out the excited levels. Taking these changes 
into account, the ground state properties are qualita- 
tively unchanged compared to Ezawa's results. We then 
present an argument showing that the degeneracy of the 
zero-energy levels is indeed conserved to all orders in per- 
turbation theory that are generated during the fRG-flow. 
The essential ingredient for this argument is the imbal- 
ance in the numbers of sublattice sites. Next we move on 
to bow-tie nanostructures which, on a bare level, also fea- 
ture zero-energy states. In this case however, the number 
of A and _B-sublattice sites is equal, and previous works 



have argued that the spins form a total singlet. Using 
the fRG we derive an effective Hamiltonian for the low- 
lying states. Now the zero-energy states are no longer 
protected and split up. The essential term in the effec- 
tive Hamiltonian that favors singlet formation is a gen- 
eralized pair-hopping term rather than a straightforward 
exchange term. Hence we conclude that only zero-energy 
states protected by sublattice number imbalance are ro- 
bust under integration of the excited states while in other 
cases, for predicting the spin ground state, the effective 
splitting needs to be compared with the interaction pa- 
rameters. 



II. EFFECTIVE ACTION 

Our aim is to describe the GNDs and related struc- 
tures within an effective theory for the low-lying single- 
particle states, in this case zero-energy-states, only. In 
this section we describe how to derive an effective the- 
ory for states near the Fermi energy by using the func- 
tional renormalization group for one-particle-irreducible 
vertices (for a derivation for bosonic field theories, see 

m)- 

We study a model described by a fermionic action 
SiW, W) of the form 

with Grassmann fcrmion fields "0 and ip depending on 
some quantum numbers (e.g. level index, Matsubara fre- 
quency, spin, etc., not written out here), the free prop- 
agator Qq containing hopping terms, chemical potential 
and Matsubara frequencies. V denotes the interaction 
which will later be assumed to be of quartic order in the 
fermion fields. We split the propagator Qo in two parts 

go = g^ + gl (s) 

and parameterise them by a matrix x 

g^ = ii-x)go (4) 
go = xgo- (5) 

At first X is arbitrary and will be specified later, e.g. 
by dividing the single-particle spectrum into low (1) and 
high energy (2) states. Then x will be a function of the 
single-particle energy e, almost zero for e smaller than 
a threshold A, and almost 1 for e > A. The partition 
function can then be split in the following form |16j 
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Z 



1 

1 

^0 



-Veff({>Al},{Vl})] 



Here, Zq = / ■0] exp [(t/;, C/f^^?/!)] is the non- 
interacting partition function, or its analogue in the case 
with superscripts. Obviously, Vcfr {V'''"}) is the ob- 

ject we are interested in: the non-trivial part of the action 
of the remaining 1-modes after the 2-inodes have been 
integrated out. Note that both types of fields, V'''" ^nd 



"0^, carry the same quantum numbers and the associa- 
tion, which degrees of freedom (e.g. high or low energy) 
they correspond to primarily is implemented through the 
choice of the cutoff x in the bare propagators. By the 
substitution ?/;^ — > ?/;^ — -!/)■*■ and -0^ —J' -0^ — -i/j-*- we get 



exp[-Vcff(Wn,{V'n)] =exp[(0;i,[c;„2]-V')] ^ / I?[Vi^0"]exp[(0^[^;o2]-V')-^({0'},W'})] 



(7) 



Now we define the effective action by 

(8) 

Here we have absorbed the integral part in ([t]) into the 
functional W"^, which under inspection turns out to be 
the generating functional for the connected Green func- 
tions with free propagator Q"^ and source-fields [f/p] '^'^ 



and [C/gJ ^ '^'^ ■ The superscript ^ indicates that this 

function includes the contribution from the 2-modes. 
Later the 2 will be replaced by an energy scale A, then 
the superscript ^ stands for "includes renormalizations 
from everything down to scale A" . 

Quite generally, the effective action derived this way 
contains arbitrarily high powers of the ^ '0-'--fields. In 
order to develop a physical picture it is most appropri- 
ate to expand the effective action in powers of the fields. 
The quadratic term then represents the renormalized free 
part, while the fourth-order term is the effective interac- 
tion [a] . Here we will not consider higher order contribu- 
tions. They are absent initially, and if the interactions 



[a] The zero-order term results only in a global shift of the energy 
and can therefore be neglected. 



are reasonably small, they should not play a decisive role. 
However, we note that this truncation issue has not been 
explored in much detail. If we now expand with 
respect to the source-fields, the quadratic part of the ef- 
fective action is given by 



(9) 



^rod is the reducible selfenergy, defined by = Gq + 
GqI^^^^Gq- From the Dyson-equation we get the relation 

^rcd = (1 - eo^=^)~^ where is the irreducible 
selfenergy. The quadratic part of the effective action is 
then 



(1 



^1 
Gox^^ 



0^; 

(10) 



In the next step we specify the matrix x- In the eigenbasis 
of C/o, X is given by 




if El < A 
if e, > A 



(11) 



with a scale-parameter A. are the eigenenergies of the 
free Hamiltonian Hq. This sharp division means that 
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degrees of freedom with > A are solely represented 
by the 2-fields, while the low-energy degrees of freedom 
with Cj < A are taken into account via the 1-fields. But in 
principle softer and also completely different definitions 
of X would be possible, resulting in different effective the- 
ories. 

In the eigenbasis of tjg the matrix M = 1 — CJoX^^ 



is not necessarily diagonal as the selfenergy can in prin- 
ciple have non-diagonal entries, e.g. in cases without 
translational invariance as the nanodiscs considered here. 
Mij — M>< would denote a component where the left in- 
dex belongs to a state i with ci > A, and the right index 
to a state j with < A. Using the cutoff-definition in 
(111, M has the structure 







> A 


< A 






M = 


> A 


1 - 


[^o]»S^> 




o]>>^>< 




(12) 




< A 






> A 




1 


< A 




M-^ = 


> A 










(13) 




L < A 





1 





We now (formally) redefine the fields by 



^4>^ = (1 -x)"^^^V^ 
-^-^^ = (1 -xV^^^-ip^ 



(14) 
(15) 



Then the quadratic part of the effective action becomes 



(16) 



with, after again using (11), 



[^efl ]» = [^'0 ^]» ' 
[5eff]<>=0 
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(17) 
(18) 
(19) 
(20) 



Next we consider the effective interactions. The quartic 
part of the effective action is given by 



(21) 



By the Dyson-series and the relation tjj'^ (fc^, fcj; /ci, ^2) 



with the two-particle- vertex 7I, it follows 



S. 



(4) _ 
off' ~ 



1 



4 

X 72 (^i^ 91. 92 



qiS 



- c 

g' [g^ 



2 \r'2 



927 



1 



^--^i^j [i + ^^g^ + ...] 



2y^2 



[1 + ^^g^ 



X7f(9i,92;9i,92) [1+g'^^' 
Here and in the rest of the section we used the Einstein 



(22) 



summation convention. Again we scale the fields by (14) 



and parameterise the propagators by the matrix x- It 
follows 



S. 



(4) 
off 



x)T [1 



X [1 
X [1 

X (1 



- S^^oX 4 
^oxS^ 4 

X)f (1 



••]/3,^ 71(91,92; 91, 92) 

[i + c;oxs^ + ... 



qiS 



(23) 



In the following we neglect the frequency-dependence 
of the selfenergy and the two-particle-vertex. We also ne- 
glect the third term in ( 20 1 and the higher orders in the 



external legs of S^'j^ which are at least linear in selfenergy 
matrix-elements that couple zero-energy-states to higher 
energy- levels. This approximation is allowed, if such ma- 
trix elements are small. In the examples described below 
this can be checked explicitly and turns out to be true. 
After this the effective action has the form 



'off ^ 



Ul^}Ves{a,P;j,5)i>^i^j (24) 
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where 



and the two-particle vertex-function 72 are 







> A 


< A 




> A 






>> 







< A 








and 



(25) 



otherwise 



(26) 



S^(fc';fc) = -Sp[^V(fc',-;fc,-)] (29) 
j^ik[, fc^; fci, fc2) = Sp [S^^Uklk',, .; fci, A:2, .)] 

-Sp [S^j^i; ■ h, fc2)[e^]^72^(fci, fc^; .)] 

-Sp [5V(fc'i,-;fci,-)eV(fc2,-;fe,-)] 

-[fci fc2] ^ [^1 ^ ^2] + [fci ^ k2,ki ^ k2] 

(30) 

in which is the full propagator and is the so called 
single-scale propagator defined by 



In (251 and ( 26 ) the zero-energy-states are fully decou- 
pled from the excited single-particle-states. Therefore we 



can map the effective action ( 24 ) into an effective Hamil- 
tonian for the zero-energy-states only. This step is pos- 
sible because we neglected the frequency dependence of 
the vertex-functions. 

The two ingredients needed for the effective action 
are the (one-particle) irreducible selfenergy and the two- 
particle-vertex. These quantities can be efficiently com- 
puted with the functional renormalization group for the 
IPI vertices [HI [18] described briefly in the next section. 



dA 



(31) 



To solve these equations we neglect the flow of the three- 
particle-vertex 7^ = and all higher vertex-functions 
and take 7^ as frequency-independent. In this approxi- 
mation is also frequency-independent. We arrive at a 
finite and closed set of flow equations that can be solved 
numerically. The truncated flow equations and the eval- 
uation of the Matsubara sums can be found in the Ap- 
pendix. 



III. FUNCTIONAL RENORMALIZATION 
GROUP 

To derive the selfenergy = Yj^ and the two-particle- 
vertexfunction 7^ = 7^ we use a functional Renormaliza- 
tion Group (fRG) scheme, with decreasing energy-cutoff 
in the free propagator. In the eigenbasis of the free 
Hamiltonian, the diagonal propagator reads 



[^0]22X^ (£2) 



[Qo\r,n (e. 



(27) 



According to the preceding section, the cutoff-function 
should be a sharp cutoff like (11), but due to numerical 



reasons we have chosen a cutoff-matrix of the form 



X^e.) = 



1 + cxp 



(/3(A-|q|) 



(28) 



where the step width of the Fermi function, 1//3, has to 
be small enough to make sure that in the end of the flow 
the zero-energy-states are not integrated out. 

The one-particle-irreducible-(lPI)-vertex-functions on 
scale A can be calculated by an infinite set of exact flow- 
equations |T71 HHJ- The equations for the selfenergy E"^ 



IV. NUMERICAL RESULTS FOR TRIGONAL 
AND BOW-TIE STRUCTURES 

Here we describe the results obtained by the numerical 
solution of the fRG flow equations for trigonal nanodiscs 
and bow-tie-shaped structures obtained by connecting 
two nanodiscs, as shown in Fig. 1. We start with the 
bare Hamiltonian as given in Eq. ([T]) for the particle- 
hole symmetric case. By integrating out the higher en- 
ergy single-particle levels down to a scale A, symmetric 
around zero energy we calculate the parameters of the ef- 
fective action for the zero-energy-states of the quadratic 
part of the bare Hamiltonian. By taking the flowing ir- 
reducible selfenergy and the IPI-vertices 7^ as fre- 
quency independent (which we already assumed in our 
truncation scheme of the vertex-functions) we can inter- 
pret them as matrix elements of an effective Hamiltonian 
for the zero-energy-sector of the free bare Hamiltonian, 



E [^0 + ^ 



a, ,a. 



I. cr, ,4i<Ti 



[72I 



»1 >»2 ''1 ''2 

c 1 ,cr2 5(^1 



(32) 



The indices ij run over all unperturbed single-particle 
states in the zero-energy sector, ct^ are the spin z- 
components. The eigenvalues of Ho + ^^ are the effective 
single-particle levels, while the second part represents 
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the effective interaction. If we assume spin-rotation- 
invariance the selfenere; yE'^ is diagonal in spin-space and 
the two-particle IPI vertex-function with a general non- 
local form can be parameterized by |18| 

72 ((a;i,cr'i), {X2,(J2)'-, {Xl,<^l), (2^2, (T2)) = 

- V {x[, X2:, Xi, X2) Sa^.a'Ja[,a2 (33) 

From the antisymmetry of 72 (fc'j^, fcj; ^ij ^2) under the 
permutations k'l ^ k'2 and fci o /c2, it follows, that the 
coupling-functions V and V obey the relation 

V [x'i,x'2]Xi,X2) = V {x'i,x'2]X2,Xi) 

^V{X'2,X'^-XI,X2). (34) 

For this reason we can simplify the effective Hamiltonian 
to 

i'l.il 
fl 

+ ^ ^it'i^«2,ii4,,^i«li,<T2°^«2.<T2a'i,<Ti (35) 

ii,if,,il,i2 

cri,cr2 

The zero-energy-states in trigonal GNDs can be cho- 
sen as eigenstates of the rotation operator -R27r/3 such 
that -R27r/3|fcj = exp(ifc)|fc,n) with k = 0, ±27r/3 and 
S\k — -|-27r/3,n) = \k — — 27r/3,n), where S is the re- 
flection operator at one symmetry axis of the nanodisc. 
Because the states \k — ±27r/3,n) are connected by a 
symmetry operation, they are degenerate in energy, while 
energy-singlets can be characterized by k=0. The num- 
ber of possible coupling functions is then reduced due 

to {iv,i2'\ii,i2) = {iv,i2'\ii,i2)5ki+k2,K,+k2, 
where the notation ij — [kj, nj) is used for the quantum 
numbers. Analogous the zero-energy-states in the bow- 
tie-shaped nanostructure are chosen as eigenstates of the 
rotation operator i?^ with i?7r|fc,n) — exp{ik)\k,n) and 
fc = 0,7r. 

First let us discuss the results for the trigonal zigzag- 
nanodiscs with iV = 2, 3 or 4 where t] is equal to the 
sublattice-imbalance Lb — La = N. Here, the first ob- 
servation by diagonalizing the quadratic part of the effec- 
tive Hamiltonian of the unperturbed zero-energy states 
is that the flow of the selfenergy remains zero for all 
zero-energy states. Hence, in these trigonal GNDs, the 
zero-energy single-particle states of the bare dispersion 
remain unsplit in the effective theory as well. For N = 2 
and = 3 particle-hole symmetry and the geometrical 
symmetry can be used to understand that there is no 
splitting. For N = 2 the two states are distinguished by 
the quantum number k — ±27r/3 and hence are degener- 
ated. Due to particle-hole symmetry they cannot move 
away from zero energy. For A'^ = 3 there is an additional 
state with A; = 0, which is also pinned to zero energy by 
particle-hole symmetry. However, for = 4 there are 



two states with fc = and one pair with k = ±27r/3. 
The two k = 0-states could be expected to split up to 
positive and negative energies, and our numerical finding 
of a robust degeneracy is surprising in the first place. In 
the next section we will give an analytical explanation 
for the protected nature of these states for arbitrary N 
and show in general under which conditions a splitting 
of the zero-energy-states will not occur. 

Besides the single-particle energies, the parameters of 
the effective interaction between the zero energy states 
have to be determined. Here, the main parts of the 
Hamiltonian ifoff can be ascribed to a direct and to an 
exchange part. These are given by 

=^ E '^n,.2«n< (36) 

il.i2 

^off - - E J^^-^-^- + (37) 

T-l .^2 
115^42 

with the matrix elements Ui^^^ = 12 ii ^^d Ji^^i^ = 

^iti2M:i2 ^^'^ spin-operators Si = lal ,.as,s'az,s' ■ In 
Fig. [5] we show the flow of Ui^^^i^ and Ji^.ij as functions 
of the flow-parameter A. As can be seen these matrix 
elements do not change drastically during the fRG-flow, 
particularly for the iV = 4 system. The initial values 
correspond to the parameters used in the analysis of 
Ezawa [H [H [H]. For N = 2, there is a simple hierar- 
chy. The largest couplings are the intraorbital repulsions 
between electrons forming a singlet in the same state. 
The next largest term is the spin exchange coupling, i.e. 
the Hund's rule coupling, and then the interorbital repul- 
sion. This hierarchy lets us already expect that for half 
filling of the zero energy levels, singly occupied orbitals 
are preferred and that the spins in this orbitals form the 
maximal total spin. For A'^ = 3 the picture remains sim- 
ilar, although now the three zero energy states consist of 
two states connected by the discrete rotational symmetry 
and another state that is strongly localized at the edges. 
Hence the couplings of this state (labeled '1') in Fig. [2] 
are larger than the ones that do not involve this narrow 
state. For a given pair of zero-energy states the hierarchy 
intraorbital repulsion > Hund's rule coupling > interor- 
bital repulsion is still visible. For A^ = 4 the picture 
is more complicated. Besides these intraorbital, Hund's 
rule and interorbital couplings there are other terms in 
HcS that show a rather mild flow as well. By this anal- 
ysis of the small nanodiscs we see that the Hartree-Fock 
analysis in [5| is already a good description of the zero- 
energy-sector. The problematic situation where a renor- 
malization or splitting of the low-energy levels competes 
with the Hund's rule interactions between these states 
does not occur. The basic aspects are captured well by 
ignoring the effects of the excited single-particle levels. 

We can go one step further and solve the effective 
Hamiltonian exactly . Written as a matrix in Fock space, 
the effective Hamiltonian H^s with all renormalized cou- 
plings included can be readily diagonalized. As a result 
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N = 2 




N = 3 




N = 4 
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FIG. 2: (color online) Flow of the coupling-functions in the 
zero-energy-sector of trigonal nanodiscs with A'^ =2,3,4 down 
to a small scale around the zero-energy-states. For the param- 
eters of the Hamiltonian |[T]) we set U = 3t and Vi — 2t. The 
kinks of the flowing coupling-functions result from integrating 
out the discrete energy- levels of Hq. In the end of the flow the 
coupling-functions can be interpretated as matrix elements of 
an effective Hamiltonian for the zero-energy-states. 



we find that in trigonal zigzag-GNDs with sizes N = 2, 
3 and 4 the ground-state spin at half fiUing is equal to 
S = Y = ^^2^"* • ^^^^ is perfectly consistent with Lieb's 
theorem. Note however that our data were obtained for 
nonzero nearest-neighbor interactions which is aheady 
outside the strict validity range of Lieb's theorem. 



0.001 




FIG. 3: (color online) Flow of the one-particle-energies (upper 
plot) and various coupling functions (lower plot) for the bow- 
tie-shaped nanostructure with rj = 2 for U = 3t, Vi — 2t. 

The situation is different in the bow-tie-shaped struc- 
tures. Here the zero-energy-states spht up in the fRG 
flow. This can be seen in the upper plot of Fig. [3] The 
splitting is rather small of the order of for our pa- 

rameters. The groundstate of the free part of the effective 
Hamiltonian is now found to be a spin-singlet, where the 
single-particle state that has been renormalized down to 
negative energy is doubly occupied. In principle, this 
degeneracy-lifting of the single-particle levels due to the 
excited levels could compete with the effective interac- 
tions between the formerly degenerate states, in particu- 
lar if there was a stronger Hund's rule coupling between 
the two single-particle levels as in the N = 2 trigonal 
nanodisc studied before. 

Hence it is instructive and important to compare the 
effective Hamiltonians of the bow-tie-shaped structure 
and the trigonal N = 2 nanodisc. We will see that the 
bow-tie ground state is indeed a singlet, but because of 
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additional interaction terms and not because of the level 
splitting in the effective free part discussed just above. In 
both cases there are two single-particle orbitals |1) and 
1 2). In the trigonal GND they can be distinguished by 
k = ±27r/3 while in the bow-tie-shape their /c-indices are 
k = 0,TT corresponding to even and odd combinations of 
the wavefunctions localized on primarily one side of the 
bow-tie whose energies are slightly split up when inte- 
grating the fRG-flow, whereat the fc = state goes down 
in energy compared to the fc = tt state. 



The two-particle spin-triplet-states have the form 
iV-t) = ^ (|1)|2) - |2)|l))xt. In trigonal N = 2 nan- 
odiscs lipt) forms the ground-state and the first excited 
states are the two degenerated singlet-states IV'sg) = 
|l)|l)Xsg and IV-sg) = |2)|2)xss- The responsible term 
in the effective Hamiltonian is the ferromagnetic Hund's 
rule- Ji2 between the zero-energy states. This J-coupling 
is also positive for the bow-tie structure, as can be seen 
in the lower plot of Fig. [3j However, the strongest 
couplings that develop during the fRG-flow are however 
pair-hopping terms denoted by J7ii22 or C/2211 of the 
form J2i j Uiijja\ ^laj^a'Cij^a in the effective Hamilto- 
nian, and the interorbital repulsions U12 — U2i- Due 
to the latter, the spin singlet states with two electrons 
in the same orbital now become compatible in energy 
with the singly occupied states in the bow-tie-shaped 
structures. But the interorbital repulsion alone is not 
sufficient to overcome the energy gain from the Hund's 
rule J12. It is the pair- hopping term that now splits the 
symmetric and antisymmetric combinations of the spin 
singlet states = ^ (IV^i,,) + 1^^%)) and jV'^p - 

^ (iV'lg) — IV'sg))) s-nd pushes the symmetric combina- 
tion IV'lg) below the Hund's rule triplet \ipt), making the 
singlet the energetically most favorable spin configura- 
tion. So, instead of being exchange-driven, the singlet 
formation is a consequence of a pair-hopping term in the 
effective Hamiltonian. Note that such a pair hopping 
term does not occur in the effective Hamiltonian of the 
trigonal nanodiscs because there it is not compatible with 
the conservation of the fc-quantum number characteriz- 
ing the rotational symmetry of the zero-energy-states. 
Furthermore, if we had done the exact diagonalization 
without taking into account the effects of the excited 
single-particle levels via the fRG, the combined effects 
of interorbital repulsion plus pair-hopping favoring the 
singlet and the intraorbital repulsion plus Hund's rule 
coupling would just equalize each other. This can be 
seen from the initial conditions for the fRG-flow of these 
couplings in Fig. [3] they all have the same values. This 
would have led us to a ground state with fluctuating to- 
tal spin. Hence the bow-tie structure is a useful example 
that including quantum corrections by higher levels into 
the low-energy Hamiltonian can make a difference. 



V. ROBUSTNESS OF THE ZERO-ENERGY 
STATES 



In the previous chapter we have seen numerically that 
in trigonal nanodiscs the zero-energy-states remain un- 
split during the fRG-flow, while the zero-energy-states in 
bow-tie-shaped structures split. This can be understood 
more generally by an analysis of the structure of the flow 
equations for the selfenergy and the two-particle-vertex, 



(38) 



7^(«'i,i2;*i>«2) =-5 E 



y^pp (*3, *4, i'i) 



13,23, 14, 

I2 («3>*4;«1:*2)72^ («i,«2;«3,«4) 
~ -^p/i(*3j i's', H, i'i) 

+ 72^ («2:*4;*l7«3)72^ (»3'«i;«4,i2) } 



(39) 



The loops on the right hand side are given in the 
appendix, Eq. (Al). Furthermore, the analysis re- 
lies on a set of properties of the electronic spec- 
trum of the free part of the Hamiltonian Hq = 

(^c] ^Ci_„ + h.c. 



Let \A) and \B) denote states that have only nonzero 
weight on A and B sublattice, respectively. Then a 
realspace-operator O shall be denoted as odd if it has 
only nonzero matrix elements of the form (A|0|i?) or 
{B\0\A). Similarly, an operator £ that has only nonzero 
matrix elements of the form or {B\£\B) is de- 

noted as even. Note, that an odd operator with real 
matrix elements is particle-hole-symmetric, i.e. invariant 
under the transformation Ui — > ^ia| with = 1 if the site 
i belongs to sublattice A and = — 1 if the site i belongs 
to sublattice B. The free Hamiltonian Hq is odd. There- 
fore the rank of the matrix Hq is at most 2La [a] and 
the nullity of Hq, which is equal to the number of zero- 
energy-states, is r] = La + Lb — rank (^Hq^ > Lb — La 

[7], which is consistent with the statements about the 
number of zero-energy-states given before. Note, that 
there are at least Lb — La zero-energy-states resulting 
from the sublattice-imbalance. Let \E) = \A) -t- \B) be 
one of the remaining 2La eigenvectors of Hq with en- 
ergy E (whereat the case = is also possible). Then 
I — i?) = |A) — |_B) is also an eigenvector of Hq with en- 
ergy —E. Therefore the remaining energies come in pairs 



[a] In the following we assume Lg > La 
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±E. The symmetric and antisymmetric linear combina- 
tions Is) = + \-E)) and \as) = - | - i?» 

are fully sublattice-polarized, i.e. |s) € A and \as) e B. 
The symmetric and antisymmetric combinations that can 
be formed from the Lb — La zero-energy-states result- 
ing from the sublattice-imbalance either vanish or are 
fully sublattice-polarized on sublattice B, because the La 
symmetric states form a basis for the states on sublattice 
A and the zero-energy-states are linearly independent on 
them. 

In the following we assume that there are no initial 
interactions of the form V{i E A,j E A;k g A, I e B) 
or V{i G B,j E B;k E B,l E A) (and cyclic), which is 
the case in the Hamiltonian Q. By this we show that 
the sclfenergy that is odd initially remains odd under the 
fRG-flow and the zero-energy degeneracy resulting from 
the sublattice-imbalance is conserved. 

The transformation matrix from the eigenbasis of Hq 
into the basis of the symmetric and antisymmetric states 
is given by 



In the following we will first show that the 



Vo 









" 


s 








as 















1 . 



(40) 



where "-f", "-" and "0" denote the sector of positive, 
negative and zero energies. Empty fields are filled with 
zero matrices of appropriate size. The transformation 
matrix from this (s/as)-space into the real space has the 
form 



(41) 







s 


as 


0" 




A 










B 









where " *" represent an arbitrary matrix block. The gen- 
eral form of odd and even matrices in the (s/as)-space 
is 





s 


as 


0" 


s 








as 




* 











* 



/ as 





s 


as 


0" 


s 








as 


















(42) 



In the eigenbasis of Hq the cutoff-matrix has the 
form 







+ 




0" 




+ 












X^- 
















(43) 



with diagonal matrices X+_|_ = X--- By transforming 
this matrix into the (s/as)-space X^/as ~ ^oX^^^ has the 
same form as £s/as Sind is therefore even. Analogously it 



follows that 
even. 



(X^)'^' and x^(x^) ^ = /3(l-x^) 



frequency-integrated single-scale propagator (see (31) 
and (A3)) is odd if we assume that the sclfen- 
ergy is odd. From the structure of the fiow equa- 
tion (38) it can be seen that the selfenergy re- 
mains odd during the fRG-flow if no two-particle- 
vertices of the form {i & A, j E A; k E A, I E B) or 
J2 {t- & B,j E B;k E B,l E A) (and cyclic) are generated 
in the flow equations ( 38 ). In a second step we show that 



these vertices will not occur, if there are no initial inter- 
actions of this form (as is assumed further above). 



Step 1: Analysis of the frequency integrated 
single-scale propagator 

If we assume that the selfenergy T,^ is odd, the ma- 
trix Hq + (x^) {x^) is odd and has a sym- 
metric spectrum. Analogue to the case of Hq we can 
transform its eigenstates into symmetrized and antisym- 
metrized states. The transformation-matrices W and V 
have the same structure as Wq and Vq- 

The first term of the frequency-integrated single-scale 
propagator ( A3 1 is 



{xY'x^ix^'y'uFU^ ixY' 
= {x^f"x^ {xT'wF,/asW^ {X^f'\, 
with the matrix F given by 



Fa,b = ( npiEa) - - ) 5aM 



(44) 



(45) 



In the eigenspace of Hq + {x^) (x^ 
F has the form 



F = 





+ 




0' 


+ 


F++ 










F__ 














the matrix 



(46) 



++ 



-F__. In the (s/as)-space F^Us = VFV^ 



and is therefore odd. 



with F. 

has the same structure as O, 
Since the remaining matrices in ( [44| are even, the first 
term of the frequency-integrated single-scale propagator 
is odd. 

The second term of the frequency-integrated single- 
scale propagator is 

{x^f'ucu^ {x^f'' 

with the matrix C given by 



Ca.l 



where ^ = (x^)'/' (x^)'''x^ (x^) 



a.b 



(47) 



1/2 . , 
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In the eigenspace of Hq + (x'^)^^^ (x^)^^^ the ma- 
trix J^"*" has the structure 





+ 







+ 




H- 


T+ 








F+ 












(48) 



with = I'l 

The matrix D ^ (x^)^^^ (x^) 
and therefore the matrix D = WDU has the structure 



•^0+ — -^0- -^+0 — 



1/2 . , 

X 



D = 





+ 







+ 


D++ 


D+- 


D+o 




D-+ 




D-o 







Do- 





(49) 



where -D++ = — -D , = —D y, 0^+ = —-Do- and 

f+o = — -D_o- So D has the same form as VOg/agV^ 
and is therefore odd. 

The structure of the matrix C in the eigenspace of 



Ho 



(X^)'^'S^ (^A)l/2 ^^^^ 







+ 







c = 


+ 


C++ 




C+0 




C-+ 













Co- 





(50) 



with C++ = -C__, C+_ = -C_+, Co+ = -Co- and 
C+o = —C-o- The matrix C^/as — VCV^ has the same 
structure as Og/as and is therefore odd. By this we see 
that the second summand of the frequency-integrated 
single-scale propagator is also odd. 

In summary it follows that the single-scale propagator 



(A3) is odd, when we assume that the selfenergy is odd. 



Step 2: Analysis of the flow equation for the 
two-particle- vertex 

In the second step of our proof, we will 
now show that no two-particle-vertices of 
the form 7^ (i G A, j G A; fc G A, Z G B) or 
^2 {i ^ B,i e B;k £ B,l £ A) (and cyclic) are 
generated in the flow equations (39). We as- 
sume that there are no initial interactions of 
the form Vii G A,j G A]k G A,l G B) or 
V{i e B,j e B;k e B,l e A) (and cyclic). From 
the structure of the flow equations ( 39 ) it can be 
seen, that no vertices of the indicated form are 



generated, if there are no Matsubara-sums (Al) 
-A — 

"pp/ph 



of the form {q e A,q' e A;k e A, k' G B) or 



I^PP/ph {q&B,q' €B-k& B, k' G A) (and cychc). 



be written in the form 

C%{q, q'- k,k')^-Y,\ Sk,k'i^^n)gt'('''r.) 



d 



dA 



J2Gtk'{^^n)Gt'('''n)- (51) 



By achieving the Matsubara-sum we get 

L J k.a L 

a.h 

AN 1/2 



u 



q.b 



111 



a,k' 
/2 



u 



,b.b' 



k,a 



6,9' 



q,b 



^a,b^a.a'Sb,b' 



b'..q' 



ix^f'w 



,1/2- 



a.'k' 



(x 



AN 1/2 



(52) 



In the (s/as)-space J^'^i^^ has the form 



/ as 





s 


as 


0' 


s 


* 






as 




















s 


as 


0" 


s 




* 




as 


* 
















I as 





s 


as 


0' 


s 








as 


















la 





s 


...as 

as 


0' 


s 








as 




















s 


,.,as 

as 


0" 


s 








as 


















j as 



,. ,as,0 





s 


as 


0" 


s 








as 


















s I as 





s 


as 


0" 






s 


as 


0' 


s 










s 








as 










as 











* 













* 





s j as 



s I as 





s 


as 


{)' 


s 








as 


















(53) 



The Matsubara-sum of the particle-hole-channel can which can be easily seen by a transformation with V. In 
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the realspace the matrix has the form 



pos 



^ pos 



,B.A 





A 


B 


A 






B 








A 


B 


A 






B 







pos 



...,A,B 



,B.B 





A 


B' 


A 






B 








A 


b' 


A 






B 







(54) 



We see that there are no expressions of the indi- 
cated form and hence no vertexfunctions of the form 
7^ (i G A, j e A; e ^, Z G B) or 
^2 {'i' ^ B,j G B;k £ B,l £ A) (and cychc) are gen- 
erated in the flow equations ( [39| . The analysis of the 
particle-particle-channel is analogous. 

In summary we have seen that when we neglect the 
frequency dependence of the vertexfunctions and trun- 
cate the fRG-equations (as usual) by setting 7^3 = 0, 
the selfenergy remains odd during the fRG-flow and 
the high zero-energy-degeneracy, following from the sub- 
lattice imbalance is conserved. This result is not re- 
stricted to a special geometry or size of the GNDs. 
For the bow-tie-shaped structure, the sublattice imbal- 
ance is zero and therefore there are no protected single- 
particle levels with zero energy. Note, that within this 
approximation the calculated selfenergy is particle-hole- 
symmetric, although the initial interactions do not have 
to be particle-hole-symmetric. A nearest-neighbor in- 
teraction like in ([T]) is for example not particle-hole- 
symmetric, but anyway the corresponding initial condi- 
tions contain no coupling-functions of the form V{i G 
A,j e A;k € AJ e B) or V{i e BJ e B;k e B,l e A) 
(and cychc). 



were obtained by neglecting the renormalizations by the 
higher levels are hence found to be valid. The protection 
of the zero-energy levels can be understood analytically 
for arbitrary N and also holds for other geometries with 
sub-lattice imbalance. If the zero-energy states occur for 
zero sublattice-imbalance, as for the bow-tie nanodisc, 
they can split up under the fRG flow, and at least in prin- 
ciple the splitting may influence the low-energy picture, 
in particular if it gets large compared to the effective 
interactions between these states. In the cases studied 
here, the splitting on the single-particle level remained 
small and was clearly dominated by interaction effects. 



In our applications of the method to trigonal nan- 
odiscs the inclusion of the excited states turned out to 
be at most a quantitative effect. Here, the large-spin 
ground state discussed in the literature is not altered 
by the renormalization. Also, the small degeneracy lift- 
ing of the single-particle levels in the bow-tie structures 
^ 10~^t ^ 2meV is in principle observable but it will be 
dominated by interaction effects. The renormalization 
of the effective interactions is also quite moderate. Nev- 
ertheless, our analysis of the bow-tie structures shows 
that integrating out the excited levels in the fRG flow 
tips the balance in the effective interactions toward a 
spin singlet rather than selecting the Hund's rule triplet. 
The fRG flow of the pair-hopping term between the ef- 
fective orbitals turns out to be stronger than that of the 
spin-exchange interaction. This nicely demonstrates the 
usefulness of the renormalization group in situations of 
competing trends. 



VI. CONCLUSIONS 

We have described a general fRG framework to to de- 
rive effective Hamiltonians for the low-lying states of 
finite-size or nanostructured lattice systems. It allows 
one to assess the influence of empty or filled single- 
particle states away from the Fermi level on the spec- 
trum and the interactions of the degrees of freedom near 
the Fermi level. With the resulting effective Hamilto- 
nian at hand, we have then performed exact diagonaliza- 
tion studies of the effective Hamiltonians in order to de- 
termine the ground-state spin of trigonal nanodiscs and 
bow-tie-shaped structures. Of course, other properties 
like transport can also be studied using the effective de- 
scription delivered by the fRG. 

The application of the fRG scheme to different smaller 
nanodiscs showed that there are two classes of nanodiscs 
(assuming nearest-neighbor hopping only) . One class has 
nonzero sublattice imbalance equal to the number of zero- 
energy states rj. Here the zero-energy single-particle lev- 
els are protected under integrating out the excited single- 
particle levels. The results in the literature[TJ UJ |B] that 



The numerical implementation of the functional renor- 
malization group scheme to these small systems is 
straightforward and one could imagine many other fields 
of applications. Another possible example with inter- 
esting low-energy states are conducting edge states of 
wires with gapped bulk spectrum, where the bulk states 
could be integrated out yielding the nontrivial effective 
description of the edge states only. However, for larger 
systems than the ones studied here, the effective inter- 
actions should be truncated in range or parameterized 
differently, otherwise the numerical effort becomes sig- 
nificant. 



We thank Sabine Andergassen, Fakher Assaad, Moto- 
hiko Ezawa, Volker Meden, Manfred Salmhofer, Jacob 
Schmiedt for useful discussions, and Bjorn Trauzettel for 
drawing our attention to graphene nanodiscs. 
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Appendix A: Evalulation of Matsubara-sums 



with 



The truncated flowequations for the ver texf unctions 
T,^ and 7^ are given in the main text, eqs. (39 1. Here 



^pp 



(Al) 



The Matsubara-sums in (38) and (Al) can be calcu- 



lated analytically. Therefore we write the Single-Scale- 
Propagator as 



X ix'y 



1/2 



ix')' 

(A2) 



Let U be the (A-dependent) transformationmatrix from 

the eigenspace of Hq + (^x^Y^^ '^^ {x^Y^^ in the re- 
alspace. The Matsubara-sum in (38) can then be trans- 



lated into a contour-integral in the complex plane and 
solved by the residue-theorem. The result is 



T.sui^^n)^p{j:[ixY'x'{ 



X npiEa) 



1/2 



u 



X") u 



[f/t {x''f'Y.^{x^)"\^{x\'u 



± ^ j v!p{Ea) if Ea = ±Eb 



(A4) 



In the contour-integral we used the function rti?(w) — 1/2 
instead of a normal fermif unction. In the Matsubara- 



sums ( Al ) we replace the free propagator Qq by the full- 



propagator (Katanin-like refinement) [19] . After the con- 
tourintegration we get 



•. 1/2 

ix^) u 



1/2 



u 



- [ix^r 
^Z>>[ix^)' 

-f^T.[ix'y 



q.a 
/2 



ix' 

1/2 



,1/2 



a,q' 



uHx^) 



u 

/2 



u 



a.b.c 



[uHx'f' 



[q ^ s] [q' O s'] 



AN 1/2 



<-'a,b.c 



(A5) 



/2 I with 



and the matrix 



1/2 



(A6) 





1 ..j.^-^/ 


nF{Et)-nF(Ea) 


n'p(Ea) 


(E^-EiY ^ 


Ea-Et 


nF(Ea)-np(Et) 


n'p(Ef,) 



if Ea 
if Ea 
if E„ 



Eb — TEc 
Eb ^ tEc 
TEc ^ Eb 



a,b 



{Ea-Et){EaTEa) ' ( E^- E^) {E^T Ea) 
^(TBc-Ba)(TBc-i5b) 



An1/2' 



if Eb = TEc^Ea 

if Ea^Eb^TEc^Ea 

(A7) 



(A3) 
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